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Abstract  In  this  paper  a  simple  crack  model  is  proposed  and  the  behavior  of  this  crack  under  stress  is 
explored.  The  crack  consists  of  an  ordinary  Griffith  crack,  but  on  either  side  of  this  crack,  at  a  distance 
tv.  exist  two  slip  planes  that  are  parallel  to  the  crack  plane.  It  is  assumed  that  slip  can  onl\  take  place 
on  these  two  planes.  Elsewhere  the  material  is  elastic.  When  it-  is  set  equal  to  zero  the  crack  becomes  a 
Bilby-Cottrell-Swinden  crack.  This  crack  model  simulates  in  a  simple  way  the  elastic  crack  tip  enclave 
model  of  a  crack.  Because  ir  has  a  finite  value  the  material  around  the  crack  tip  is  elastic.  The  crack  is 
considered  to  be  stressed  in  either  mode  II  (plane  strain  shear)  or  mode  111  (anti-plane  strain  shear).  It 
is  found  that  for  a  virgin,  stationary  crack  the  stress  intensity  factor  at  the  crack  tip  is  equal  to  the 
conventional  stress  intensity  factor  when  the  stress  is  raised  under  a  monotonically  increasing  load. 
However,  when  the  crack  tip  advances  under  such  a  load  the  crack  tip  stress  intensity  factor  is  smaller 
than  the  conventional  stress  intensity  factor.  The  fracture  stress  is  proportional  to  the  surface  energy  of 
the  solid  raised  to  a  power.  In  general,  this  power  is  not  equal  to  one  half.  For  cyclic  loading  by  using 
qualitative  arguments  it  is  shown  that  the  crack  can  grow  an  incremental  distance  each  cycle,  and  the 
growth  law  is  a  fourth  power  Paris  equation. 

Resume— Dans  cet  article,  nous  proposons  un  modele  simple  pour  une  fissure  dont  nous  etudions  le 
comportement  sous  contrainte.  Cette  fissure  est  une  fissure  ordinaire  de  Griffith,  mais  il  existe  de  chaque 
cote  de  la  fissure,  a  une  distance  tv.  deux  plans  de  glissement  parallels  au  plan  de  la  fissure.  Nous 
supposons  que  le  glissement  ne  peut  se  produire  que  sur  ces  deux  plans.  Partoul  aillcurs.  Ic  material)  est 
elastique.  Lorsque  ir  egale  zero,  la  fissure  devient  une  fissure  de  Bilby.  Cottrell  et  Swinden.  Ce  modele 
de  fissure  permet  de  simuler  simplement  le  modele  d'enclave  d'une  extremite  de  fissure  elastique.  ic  ayant 
une  valeur  finie,  le  materiau  autour  de  l'extremite  de  la  fissure  est  elastique.  Nous  considerons  que  la  fissure 
est  mise  sous  contrainte  selon  le  mode  II  (cisaillement  plan)  ou  le  mode  III  (cisaillement  anti-plan).  Dans 
le  cas  d'une  fissure  vierge  stationnaire.  le  facteur  d'intensite  de  la  contrainte  a  l'extremite  de  la  fissure  est 
egal  au  facteur  d  inlensite  de  contrainte  classique  lorsque  Ton  augmente  la  contrainte  de  maniere 
monotone.  Cependant,  lorsque  l'extremite  de  la  fissure  avance  sous  1'effet  d'une  telle  contrainte.  le  facteur 
d’intensite  de  la  contrainte  a  l'extremite  de  la  fissure  est  plus  petit  que  le  facteur  d'intensite  de  contrainte 
classique.  La  contrainte  de  rupture  est  proportionnelle  a  une  puissance  de  1'energie  superlicielle  du  solide: 
en  general,  cette  puissance  n’est  pas  egale  a  un  demi.  Dans  le  cas  d'une  charge  cyclique.  dcs  arguments 
qualitatifs  montrent  que  la  fissure  peut  se  propager  sur  une  distance  qui  augmente  a  chaque  cycle,  et  que 
la  loi  de  croissance  est  une  equation  de  Paris  de  degre  quatre. 

Zusammenfassung — Ein  einfaches  RiBmodell  wird  vorgeschlagen,  Insbesondere  wird  das  Vcrhalten  des 
Risses  unter  einer  angelegten  Spannung  untersucht.  Der  RiB  ist  ein  gewohnlicher  Griffith-RiB.  Auf  beiden 
Seiten  dieses  Risses  jedoch  finden  sich  zwei  Gleitebenen  parallel  zur  RiBebene  im  Abstand  »\  Es  wird 
angenommen,  daB  Gleitung  nur  auf  diesen  beiden  Ebenen  ablauft,  ansonsten  ist  das  Material  elastisch. 
Mit  tv  =0  erhalt  man  einen  RiB  entsprechend  Bilby-Cottrell-Swinden.  Das  vorgelegt  Modell  beschreibt 
auf  einfache  Weise  das  Enklavenmodell  einer  elastischen  RiBspitze.  Da  tv  einen  endlicher.  Wert  besitzt, 
ist  die  RiBspitze  elastische.  Es  wird  angenommen,  daB  der  RiB  entweder  in  Mode  II  (ebenc  Scherung)  oder 
Mode  III  (anti-ebene  Scherung)  SpannungsmaBig  belastet  wird.  Bei  einem  frischen  stationaren  RiB  gleicht 
der  Spannungsintensitatsfaktor  an  der  RiBspitze  dem  konventionellen  Spannungsintensitatsfaktor.  wenn 
die  Spannung  mit  einer  monoton  ansteigenden  Last  erhoht  wird.  Wenn  sich  jedoch  die  RiBspitze  unter 
einer  solchen  Last  ausbreitet,  dann  wird  der  Spannungsintensitatsfaktor  kleiner  als  der  konventionelle. 
Die  Bruchspannung  ist  proportional  zur  Oberflachenenergie  des  Festkorpers  hoch  einer  Potenz,  die  nicht 
ein  halb  ist.  Qualitative  Argumente  zeigen,  daB  der  RiB  unter  zyklischer  Belastung  bei  jedem  Zyklus  ein 
wenig  wachsen  kann:  das  Wachstum  wird  durch  eine  Parisgleichung  der  vierten  Potenz  beschrieben. 


INTRODUCTION 

In  this  paper  a  very  simple  crack  model  is  considered 
in  order  to  elicit  the  fundamental  features  of  a  crack 
interacting  with  a  deformation  field  generated  by  the 

tSupported  in  part  by  U.S.  Army  Research  Office. 


stress  concentration  at  the  crack  lip,  and  to  confirm 
results  of  an  earlier  more  complex  elastic  enclave 
model  [1,2].  In  our  model,  we  assume  that  slip  can 
occur  on  two  slip  planes,  parallel  to  the  cleavage 
plane,  but  displaced  from  it  by  distance  ir  in  each 
direction.  Figs  1  and  2.  Our  work  will  be  couched  in 
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_ Slip  Plane _ 

— cr°cfc _ i—  _ _ , 

Plon*  w 

_ i _ _ 

Sl<p  Plane 
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Fig.  I  Double  slip  plane  model.  The  crack  is  assumed  to  be¬ 
an  atomically  sharp  slip  crack  on  the  negative  v-axis  with 
the  tip  at  the  origin  Two  planes  parallel  to  the  crack  plane 
are  permitted  for  dislocation  generation  and  slip.  The  planes 
are  assumed  to  be  a  distance  n  front  the  crack  plane.  Mode- 
Ill  or  II  will  be-  assumed  throughout  for  purposes  of 
analysis. 

terms  of  the  shielding  afforded  a  crack  by  the  defor¬ 
mation  field  generated  by  the  stress  concentration  at 
the  crack  tip. 

For  our  simple  model,  we  will  show  that  (I)  the 
initial  response  of  a  crack  to  the  generation  of  a 
deformation  field  on  the  two  slip  planes  is  essentially 
nonshielding  in  character;  (2)  that  as  the  crack  prop¬ 
agates,  the  stress  intensity  at  the  crack  tip.  A,  is 
shielded  from  the  external  stress  intensity.  A',  that  is 
A'  >  A',;  and  (3)  that  cyclic  loading  with  work  hard¬ 
ening  can  lead  to  crack  growth. 

These  results  are  based  upon  a  particular  form  of 
one  dimensional  slip,  and  we  can  expect  the  results 
will  be  modified  if  slip  occurs  on  planes  which 
intersect  the  cleavage  plane.  We  will  discuss  these 
modifications  in  a  qualitative  way,  but  the  simple 
model  nevertheless  leads  to  results  which  highlight 
how  the  behavior  of  materials  with  real  dislocations 
will  differ  from  that  of  a  plastic  continuum.  We  also 
note  here  that  our  one  dimensional  model  has  certain 
similarities  to  the  early  BCS  model  [3]  of  fracture,  but 
the  main  physical  picture  resulting  in  our  case  is 
different  in  important  ways  from  the  BCS  results.  We 
will  work  in  mode  III  and  mode  II  because  the 
mathematics  will  be  accessible  to  an  analytic  ap¬ 
proach  (especially  in  mode  III),  but  the  general 
features  are  expected  to  apply  to  the  more  important 
mode  I,  as  well.  Finally,  we  note  that  our  paper 
addresses  the  behavior  of  a  crack  whose  tip  is  atom¬ 
ically  sharp. 

THE  CRACK  TIP  STRESS  INTENSITY  FACTOR 

In  an  earlier  paper  [4],  by  one  of  us,  we  showed  for 
a  set  of  dislocations  in  mode  III  surrounding  a  sharp 
crack,  that 

G  =  K,  +  X  Rj  ~  ( I ) 

d  2/i 

where  g  is  the  force  per  unit  length  exerted  on  the 
defect,  here  referred  to  as  t  for  the  crack,  d  for  a 
dislocation,  and  G  is  the  force  on  the  conglomerate 
composed  of  crack  plus  dislocation  distribution.  G 
and  g  are  complex  functions  of  v  +  iy,  the  crack  is 
supposed  to  lie  along  the  negative  .v-axis  with  the  tip 
at  the  origin,  the  real  part  of  g  or  G  is  the  v- 
component  of  the  force,  and  the  imaginary  part  the 


Slip  Zone 

i 

ciiiii 


T  TTT 

x  *  0 

Fig.  2.  Double  slip  plane  model.  Dislocations  are  distributed 
symmetrically  on  Ihe  two  planes. 

i  -component,  g  is  the  complex  conjugate  of  g.  fi  is  the 
shear  modulus.  The  physical  interpretation  of  equa¬ 
tion  ( I )  is  that  the  total  force  G  =  A’  2  g  in  terms  of 
the  far  field  external  stress  intensity  factor  (small  scale 
yielding)  is  linearly  composed  of  the  separate  forces 
on  crack  and  dislocations.  Equation  (1)  is  valid  only 
when  the  dislocation  distribution  is  symmetric  (the 
situation  considered  in  this  paper)  about  the  v-axis. 
The  local  force  on  the  crack  by  itself  is  given  by 


.?,=  A’;  2  ft 
A,  =  A  -  A',, 


s',  is  the  position  of  the  /th  dislocation  in  the  complex 
plane.  A',,  is  the  screening  contribution  to  A'  at  the 
crack  tip  due  to  the  dislocation  distribution.  The 
physical  meaning  of  equation  (2)  is  that  the  crack  is 
either  screened  (or  antiscreened)  from  the  externally 
imposed  stress  represented  by  the  external  A-field, 
depending  on  the  sign  of  />,. 

In  Fig.  2,  since  the  .v -component  of  force  on  a 
dislocation  is  given  by  Re(g(/)  =  x where  x  is  the 
shear  stress  on  the  slip  plane,  we  can  also  write 

&  +  ITA=&  +  T(.v)fl(.v)d.v  =  A:/2g.  (3) 

Flere  we  have  replaced  the  discrete  dislocations  by  a 
continuous  linear  distribution  df>,  =  B(.x)  dv.  The  in- 


Fig.  3.  Path  for  /-integral.  T  is  a  path  outside  the  crack  and 
all  dislocations,  y  encloses  the  crack  tip,  and  ,Sj  and  ,S\ 
enclose  the  two  symmetric  dislocation  distributions  on  the 
parallel  slip  planes. 
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tegration  is  carried  out  on  only  one  slip  plane,  but 
in  Eig.  (2),  we  assume  that  the  distribution  on  each 
slip  plane  is  identical,  accounting  for  the  factor  2 
before  the  integration  in  equation  (2). 

Equation  (2)  is  easily  generalized  to  mode  II  and 
mode  I  by  taking  the  ./-integral  (the  ./-integral  of  Rice 
[5]  which  is  identical  in  the  elastic  case  to  the  /■', 
integral  proposed  earlier  by  Eshelby  [6])  path  shown 
in  Eig.  2  (This  is  in  fact  what  was  done  in  Ref.  [4].) 
The  external  path  enclosing  crack  and  dislocations 
yields  the  total  crack  extension  force.  G  =  K'2  JE  /i, 
where  .K  is  a  mode  dependent  factor  given  below. 
is  K]  2.X ft.  and  the  path  around  the  dislocation 
distribution  is  given  by 

c*u  \ 

IV dr  —  f—  d.v  ).  (4) 

‘  v  1 

Thus  we  can  replace  (2)  by  the  more  general 


r(.v)B(.v)d.v  =  K:  2  Jf/i 


fl  Modelll 
ll-v  Mode  II. 


(5) 


It  should  be  noted  in  equation  (2)  and  (5)  that  we 
are  considering  only  the  case  when  a  continuous 
dislocation  distribution  fl(.v)  does  not  have  a  stress 
singularity  associated  with  it.  This  is  the  case  when 
r(.v)  is  tinite  everywhere  on  the  slip  planes.  If  stress 
singularities  exist  at  an  end  of  a  dislocation  distribu¬ 
tion,  as  they  would  if  dislocations  were  piled  up 
against  an  obstacle,  the  integrals  of  (3)  and  (5)  must 
be  evaluated  with  great  care.  For  example,  in  that 
case,  a  finite  contribution  is  made  to  the  7-integral 
from  the  tips  of  the  distribution,  whereas  when  r(.v) 
is  everywhere  finite,  the  contribution  from  the  tips  is 
zero  [7], 

It  should  be  noted  that  equation  (2)  and  (5)  can  be 
rewritten  in  terms  of  the  (total)  dislocation  produced 
shielding  or  antishielding  stress  intensity  factor,  L,  of 
Rice  and  Thomson  [8).  Equation  (2)  then  becomes 


dislocation  disti ibution  on  the  discrete  slip  planes  of 
Figs  1  and  2.  but  require  the  crack  to  remain  station¬ 
ary.  Specifically.  (I)  we  place  a  slit  crack  in  an 
otherwise  perfect  elastic  medium:  (2l  we  allow 
sources  to  operate  anywhere  on  each  of  the  two  slip 
planes  so  that  identical  distributions  are  produced  on 
each.  It  is  a  property  of  these  distributions  that  equal 
numbers  of  dislocations  of  positive  and  negative 
burgers  vector  are  produced  on  each  slip  plane 
because  a  dislocation  source  always  produces  dis¬ 
locations  in  pairs.  Thus 

B(v)d.v  =(): 

(3)  we  allow  the  dislocation  distribution  to  come  into 
force  equilibrium,  so  that  the  total  force  acting  on  a 
dislocation  is  zero.  The  elastic  portion  of  the  force 
will  be  made  up  of  the  total  elastic  shear  stress,  r, 
acting  at  the  dislocation  in  question  from  all  the  other 
dislocations  and  the  stress  from  the  crack.  There  will 
also  be  a  force  due  to  lattice  friction,  due  to  the 
effective  lattice  friction  stress  which  we  call  a,.  Thus 
we  can  write  in  equilibrium  that 

t(.v,)=  -<r,(.v,).  (10) 

If  it,  is  a  simple  constant,  then  the  crack  model  would 
simulate  a  perfectly  plastic  solid.  If  a,  increases  with 
the  displacement.  D,  across  a  slip  plane  where 

/>(v)  =  j  fl(.v)d.v  (II) 

then  the  crack  model  would  simulate  work  hard¬ 
ening.  (The  reader  will  note  that  true  “work  hard¬ 
ening"  is  not  really  consistent  with  our  model,  be¬ 
cause  the  interaction  with  the  other  dislocations  is 
presumably  included  in  equation  (3)  or  (5).  In  true 
work  hardening,  "redundant"  dislocations  of  op¬ 
posite  sign  are  all  important.  The  net  Burgers  vector 
of  these  redundant  dislocations  in  any  small  volume 
element  is  equal  to  zero,  and  consequently  they  do 


K,  =  K  +  L  (6) 

where  L  has  a  positive  value  in  the  antishielding  case 
and  a  negative  value  in  the  shielding  situation.  In 
equation  (5)  let  the  integral  term  be  set  equal  to 

I  r(.v  )B(.v)  d.v  =  j—  /  (7) 

where  /  is  defined  by  this  equation.  Equation  (5)  can 
be  written  as 

K,=  K(\-I/K')'2.  (8) 

Combining  equations  (6)  and  (8)  gives 

L  =  -K+(K2 -l)'\  (9) 

SIaMoNaRA  CRACK 

We  now  explore  the  result  of  a  sequence  of  events 
in  which  we  allow  a  sharp  slip  crack  to  generate  a 


not  appear  in  a  dislocation  distribution  function 
B(x).  We  therefore  have  to  assume  for  this  purpose 
that  there  are  other  effects  which  increase  <r,.  perhaps 
lattice  debris,  other  dislocations  in  walls,  etc.,  which 
we  do  not  otherwise  explicitly  countenance.  (How¬ 
ever,  making  a,  depend  upon  D  does  give  us  an 
analytic  way  to  bring  the  real  effects  of  deformation 
into  our  simplified  model.) 

We  now  calculate  the  total  elastic  force  on  the 
dislocations  I'fum  equation  (3)  or  (5).  This  is  given  by 


r(.v)fl(.v)d.v 


<r,B(.v)d.v. 


(12) 


In  the  spirit  of  the  last  paragraph,  we  now  assume 
that  a,  =  (T,(Z)).  Since  B(.v)d.v  =d D,  we  then  have 


(T,(.v)fi(.v)d.y  = 


<!,(D)dD=Q.  (13) 
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Likewise,  /  given  in  equation  (7)  is  zero.  Substituting 
for  Igj  in  (3)  or  (5),  we  then  have 


A;  A'2 
2  JTp  ~  2  Jfp 


(14) 


(or  from  equations  (8)  and  (9)  A,  =  A  and  L  =  0),  amt 
the  crack  tip  is  not  screened  by  the  deformation  zone. 
Although  this  result  is  derived  for  two  parallel  slip 
planes,  it  is  also  valid  for  an  arbitrary  number  of  slip 
planes  parallel  to  the  cleavage  plane  in  which  sym¬ 
metry  about  the  v  axis  is  preserved. 

This  result  is  so  striking  that  it  deserves  additional 
discussion  and  verification.  (This  result,  as  pointed 
out  by  a  reviewer,  is  obvious  from  the  original 
interpretation  of  Eshelby’s  [6]  A,  integral.  A,  measures 
the  resultant  on  the  crack  tip  of  the  forces  due  to  the 
elastic  field  on  all  the  dislocations,  which  is  opposite 
to  the  traction  stress.  With  equal  numbers  of  dis¬ 
locations  of  opposite  sign,  and  each  opposed  by  the 
traction  stress,  only  the  crack  extension  force  re¬ 
mains.)  Suppose  we  have  only  one  pair  of  dis¬ 
locations  of  opposite  sign  on  each  of  the  slip  planes 
of  Fig.  4,  and  assume  mode  III.  Then  we  explore  the 
possibility  of  setting  the  dislocation  contribution  to  A 
equal  to  zero 


0  =  A',  =  pb  Ret 


pb 


.  _(  cos^  y  sinO,  -  cos  ^  N/sin  0, 
2iru\  2  2 


(15) 


For  any  value  of  0,  (where  we  place  />,)  in  the  range 
0  <  0,  <  7t,  there  is  another  and  different  value  for  0: 
(where  we  place  b2)  for  which  Kd  =  0  (except  for  the 
single  angle  where  only  one  value  of  0  satisfies  ( 1 5), 
6  =  2  tan"1  l/y/i).  From  equations  (3)  and  (4)  when 
A,  =  0 


Re(gl+g2)  =  0.  (16) 

and  if  there  exists  an  effective  lattice  friction  stress,  of, 
which  is  the  same  for  each  dislocation,  since 
b,  =  —b2,  we  have  constructed  an  equilibrium 
configuration  for  which  there  is  no  screening.  Equa¬ 
tion  (14)  is  thus  confirmed.  We  note  that  as  tv-»0  in 
equation  (15)  the  argument  degenerates. 

The  nonshielding  result  of  equation  (14)  is  a  special 
case  for  slip  planes  which  are  parallel  to  the  cleavage 
plane.  However,  even  in  cases  where  this  geometry  is 
not  approximately  a  valid  picture,  we  retain  some 
features  of  the  nonshielding  result.  The  reason  is  that 


w 


1 _ 
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Fig  4  Schematic  figure  for  two  pairs  of  dislocations  in 
non-shielding  configuration.  Coordinates  of  dislocations  are 
labelled. 


Fig.  5.  Friction  law.  Dashed  lines  show  the  law  obeyed  lor 
reversed  slip  following  forward  slip  with  work  hardening. 
The  full  line  for  negative  slip  indicates  a  fully  reversible  slip 
law. 


in  all  cases,  dislocations  will  initially  be  formed  in 
pairs,  and  the  pair  dipoles  will  be  so  oriented  that  the 
antiscreening  dislocation  will  be  closest  to  the  crack, 
and  will  be  attracted  toward  it.  Hence,  the  initial 
effect  of  loading  a  stationary  crack  will  not  be  to 
shield  it.  but  shielding  can  only  occur  as  a  later  stage 
of  the  total  process.  We  turn  now  to  a  discussion  of 
possible  later  stages. 


GROWING  CRACK 

Let  us  suppose  we  have  a  stationary  crack  with  two 
active  double  slip  planes  as  described  in  the  previous 
section  in  which  Kd  =  0  and  A  =  A,,  and  that  we 
begin  to  increase  A  in  adiabatic  fashion.  When  the 
value  of  A,  reaches  the  Griffith  critical  value  for  crack 
propagation  in  a  perfectly  brittle  solid,  the  crack 
propagates.  However,  when  this  happens,  the  re¬ 
lations  between  the  elastic  stresses,  r  and  a ,  which  are 
responsible  for  establishing  the  equilibrium  of  the 
dislocation  configuration  become  altered,  and  the 
nonshielding  solution,  equation  (10),  breaks  up.  In 
effect,  the  antishielding  dislocations  undergo  reversed 
slip  if  they  are  carried  along  with  the  propagating 
crack,  and  af  for  these  dislocations  then  changes  sign. 
See  Fig.  (5).  This  reversal  of  o,  means  that  the  crack 
will  leave  the  antishielding  component  behind,  and 
the  shielding  dislocations  will  be  pushed  ahead  of  it. 
The  net  effect  is  that  only  the  shielding  dislocations 
are  left  in  the  distribution.  The  sum  in  equation  (3) 
to  obtain  Kd  is  then  finite  and  positive,  and  A  >  A,  so 
that  the  crack  tip  is  shielded. 

That  antishielding  dislocations  cannot  be  carried 
along  by  a  growing  crack  can  be  understood  from  the 
following  considerations.  If  antishielding  dislocations 
were  carried  along  and  thus  move  in  the  same 
direction  as  the  shielding  dislocations,  the  stress  on 
the  slip  plane  must  change  discontinuously  from  the 
value  <r,  to  the  value  —  a,  if  there  is  no  gap  on  the  slip 
plane  between  the  shielding  and  the  antishielding 
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Fig.  6.  Dislocation  distribution,  (a)  Schematic  dislocations 
for  one  of  a  pair  of  parallel  slip  planes.  The  distribution  near 
.v  =  h'  falls  to  zero  smoothly,  (b)  For  the  purposes  of 
analyses,  we  assume  the  dislocation  distribution  is  on  the 
crack  plane,  and  that  the  closest  approach  of  the  dislocation 
distribution  to  the  crack  tip  is  ic. 

dislocations.  But  a  stress  discontinuity  of  this  type 
can  be  produced  only  by  a  logarithmically  divergent 
dislocation  distribution  which  has  the  same  sign  on 
either  side  of  the  discontinuity.  (This  can  be  shown 
using  the  usual  Hilbert  transform  integral  equations 
that  appear  in,  for  example,  the  BCS  crack  theory.) 
Hence  no  solution  exists  if  no  gap  exists.  But  suppose 
a  gap  does  exist.  The  antishielding  dislocation  must 
trail  after  the  shielding  dislocations.  If  they  were  in 
the  forward  position  they  would  be  pulled  into  the 
shielding  dislocations  because  of  the  attracti  le  stress 
fields  of  both  the  crack  tip  and  the  shielding  dis¬ 
locations.  Because  the  friction  stress  is  supposed  to 
have  a  reversed  sign  in  the  region  of  the  antishielding 
dislocations,  the.e  dislocations,  too,  according  to 
equation  (3)  or  (7),  are  shielding  dislocations  and 
cause  a  reduction  in  the  value  of  K,.  But  the  anti¬ 
shielding  stress  intensity  factor  L  of  the  antishielding 
dislocations  cannot  change  their  sign  according  to  the 
equations  of  Rice  and  Thomson  [5].  Thus  our  sup¬ 
position  that  there  is  a  zone  of  antishielding  dis¬ 
locations,  with  a  reversed  friction  stress  that  is  carried 
along  with  the  crack  tip  through  the  attraction  of 
these  dislocations  towards  the  shielding  dislocations, 
has  led  to  an  inconsistent  result  and  the  assumption 
must  be  incorrect. 

We  will  now  carry  out  an  analysis  of  the  shielding 
to  show  how  it  depends  in  detail  on  the  material 
parameters.  We  shall  assume  the  distance,  w,  is  small 
compared  to  the  length  over  which  the  dislocations 
are  distributed.  In  this  case,  the  closest  approach  of 
any  dislocation  to  the  crack  tip  will  be  tv.  Further, 
except  for  those  dislocations  closest  to  the  tip,  the 
angular  part  of  the  crack  is  negligible,  and  0  can  be 
taken  0=0  to  good  approximation.  Thus  we  can 
approximate  the  entire  distribution  and  its  effect  on 
the  crack  by  replacing  the  pair  of  slip  planes  with  a 
single  line  of  dislocations  on  the  cleavage  plane,  but 
where  the  distribution  begins  at  .v  =  tv.  That  is,  we 


assume  a  "dislocation  free  zone"  in  front  of  the 
crack,  to  correspond  to  the  slip  plane  separation,  u  . 
See  Fig.  6a,  b.  This  strategem  allows  us  to  carry  out 
a  rigorous  analysis. 

From  Ref.  (4)  the  force  on  an  arbitrary  dislocation 
on  the  slip  line  is 


Kj  = 


Kb 

V  2  7I.Y 


4ti.v 


+  1 


Hb  b  fx'  I 

V  x 


2  71 


v  —  v 


(17) 


.v  is  the  distance  on  the  cleavage  plane  from  the  crack 
tip  to  the  dislocation,  and  a  sum  is  taken  over  all 
other  dislocations  at  .v'  of  Burgers  vector  />'.  Adop¬ 
ting  a  continuum  limit,  we  set  db  =  B(x)dx.  and 
write 


I  dg  K  l*B(x)  d.v 

B  d  v  v  2nx  4tt.v 


V  B(.x  ') 


-Of(x) 


(18) 


ff/(.v)  on  the  right  is  the  lattice  friction  stress,  and  with 
it  (18)  is  the  condition  for  equilibrium  to  be  satisfied. 
We  shall  take  <rr  as  a  constant  in  the  following, 
because  otherwise  the  mathematics  becomes  intrac¬ 
table.  The  second  term  is  of  lower  order  than  the 
others,  and  must  be  dropped.  We  note  that  this 
dropout  of  the  self-image  term  demonstrates  a  phys¬ 
ical  difference  between  continuum  theories  and  dis¬ 
crete  theories.  The  self-image  term  is  known  to  be 
dominant  at  distances  of  order  b  from  the  crack  tip, 
and  thus  all  continuum  theories,  because  they  lose 
this  scale  parameter,  are  essentially  helpless  to  discuss 
atomic  phenomena  at  the  crack  tip.  Equation  (18) 
can  be  put  in  the  canonical  integral  equation  form 


P(x’)  =  ^  B(x').  (19) 

This  is  the  standard  singular  integral  discussed  by 
Muskhelishvili  [9]  and  we  refer  the  reader  to  the 
paper  of  Head  and  Louat  [10]  for  its  solution.  This 
problem  has  been  solved  recently  in  a  different  way 
by  Chang  and  Ohr  [11],  After  carrying  out  the 
analysis  described  below,  we  discovered  that 
Majumdar  and  Burns  [12]  in  a  paper  shortly  to  be 
published  have  carried  out  the  solution  by  a  tech¬ 
nique  very  close  to  our  own.  However,  since  we  are 
interested  in  the  analytic  results  which  can  be  ob¬ 
tained  on  the  assumption  that  »v  R,  we  shall  indi¬ 
cate  our  steps  in  summary  form  below. 

The  uniqueness  relation  of  Head  and  Louat  [10]  for 
a  distribution  B(x)  which  has  zero  values  at  it-  and  R 


47s 


WI  I.RTMAN  vi  ul.:  DOUBIF  Sl.ll>  PI.ANK  CRACK  MODI  I 


is  given  by 

A  =  -  N  2n  it,  /,  /, 


'  (  v  -  n  )(/?  -v) 


d.v  '  -  v  2  R 


+  — —  In  4  R  h- 
■v  * 


,  r  dv 

/ ;  =  — =  71. 

X,  (  V  -  w){R  -  ,v  ) 


(20) 


The  integral  equation  is  inverted  by  the  equation 


2  /( v  -  ii )( /?  -  a  )  fA7,(.v) 

fl(A)=-—  / - - V — =-  +  rr;/4 

V  v  (  v  2n 

,  f*  I  d.v' 

/,  =  -^=== - ;  =  0 

J»  v  (V  -\v)(R  a  )  v  -  v 

/=r  /  ?  dv- 

Jh  \]{.\'-\y)(R  -  v'j.v -a' 


(21) 


Using  the  uniqueness  relation.  (20),  and  carrying  out 
the  integration  to  find  A',  from  (3),  we  have 


A,  =  A  - 


v  2rt  J 


'  fl(.v )  d.v 


(22) 


Combining  the  previous  equations  we  find 


*=1+i^r 

K  n-h  , 


A(v) 


V'(A  -»')(/?  -«  )^d.v.  (23) 


We  carry  out  the  integrations  indicated  in  the  limit 
R  P  »•  by  elementary  integrals  and  find  to  lowest 
order  in  wIR 


(In  4  R/w  +  4/3) 

_3/2  _ 

A,  =  — ~r~  ai  v' M'  x  (*n  4  R/ »  4-  4/3).  (24) 
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This  result  is  equivalent  to  that  obtained  by  Ma- 
jumdar  and  Burns  [12].  Its  form  in  equation  (24)  is 
particularly  useful  to  us  because  of  its  simple  analytic 
character. 

It  should  be  noted  that  the  result  given  by  equation 
(24)  says  that  the  total  integrated  burgers  vector  of 
the  dislocations  ahead  of  the  crack  tip  is  the  maxi¬ 
mum  possible  subject  to  the  condition  that  at  a 
distance  w  ahead  of  the  crack  tip  the  stress  is  equal 
to  Of.  It  is,  of  course,  possible  to  find  a  solution  in 
which  a  smaller  number  of  dislocations  are  present 
ahead  of  the  crack  tip.  The  value  of  K,  then  is  larger 
than  that  given  by  equation  (24).  Bu.  in  this  situation 
the  stress  would  be  larger  than  o,  at  the  distance  n. 
If  more  dislocations  were  present,  the  dislocation 
distribution  would  either  have  to  start  at  a  distance 
smaller  than  »■,  or  if  the  dislocations  were  restricted 
to  start  at  w,  an  infinite  stress  at  this  position  is 
required  to  prevent  the  dislocations  from  ap¬ 
proaching  closer  to  the  crack  tip. 


K,  _  3  Iw 
~K~2n\jR 


We  note  that  equation  (24)  contains  the  shielding 
physics  which  the  BUS  model  does  not.  f  or  example, 
if  we  let  w  »()  as  in  the  BCS  model.  A  ►(),  wherein 
we  lose  the  A  field  of  the  underlying  crack  entirely, 
l  or  us.  u’  will  retain  the  meaning  of  the  separation 
between  slip  planes.  We  note  that  for  slip  planes 
actually  separated  from  the  cleavage  plane,  that  is,  as 
shown  in  Fig.  1  rather  than  Fig,  6b.  the  solution  very 
near  the  tip  will  be  somewhat  dillerent  from  the 
cleavage  plane  model  of  equation  (23).  The  distribu¬ 
tion  near  the  crack  tip  will  not  cut  off  as  sharply  as 
in  (21).  but  this  dilference  is  expected  to  be  a  minor 
quantitative  effect  rather  than  a  qualitative  one. 

To  incorporate  work  hardening  directly  into  the 
dislocation  equilibrium  equations.  (IS),  o,  would 
have  to  be  a  function  of  the  local  density  fl(.v),  and 
there  would  be  no  analytic  solution  of  the  integral 
equation.  Flowever,  if  we  allow  a,  to  be  a  function  of 
the  total  dislocation  content  of  the  slip  plane,  we  can 
then  incorporate  a  simple  work  hardening  law  into 
our  results.  Thus,  we  shall  set 


«t  =  o„(D*  D,r 
(a  )  d.v 


*R 

0*  =  B(v 

4/  M 


(25) 


n„  and  D„  are  constants,  and  m  is  the  work  hardening 
exponent.  0  <  m  <  I. 

To  the  lowest  order  in  ic  R.  for  mode  III.  we 
integrate  fi(.v)  to  obtain 


(26) 


Then,  combining  equation  (24).  (25)  and  (26)  we 
obtain 


•  «t)  Am 

li 


_ _ \  2/)„/i  _ 

o,\  ’"'(In  4  A  'ic  +4  3)' . . 


(27) 


From  equation  (27)  and  our  previous  discussion,  we 
can  now  predict  the  overall  behavior  of  a  bare  crack 
when  at  time  /  =  0.  we  load  it  up  and  allow  defor¬ 
mation  to  proceed  in  the  vicinity  of  the  crack  tip.  We 
must  first  assume  that  the  crack  will  propagate 
whenever  A,  exceeds  a  critical  value  for  cleavage  of  a 
pure  brittle  crack  in  the  solid.  A,,,.  In  an  uncom¬ 
plicated  situation,  we  must  reasonably  assume  this 
critical  value  is  the  Griffith  value 


Ki.  =  v  4-//i  .*  (28) 

where  y  is  the  intrinsic  surface  energy  of  the  solid  At 
time  I  =  0,  the  crack  is  loaded,  and  dislocation  pairs 
are  created  on  the  slip  planes.  For  a  short  transient 
period,  the  pairs  remain  in  the  vicinity  of  the  crack 


WFF.RTMAN  cl  al\.  DOUBLE  SLIP  PLAN!  C  RAC  K  MODEL 


tip.  and  the  A',  remains  near  zero.  The  external  A  is 
transmitted  directly  to  the  crack  tip  unshielded.  The 
crack  will  accelerate  as  soon  as  A,.  is  achieved,  the 
shielding  dislocations  will  be  shoved  ahead  of  the 
crack,  and  antishielding  ones  will  be  sloughed  olT  at 
the  rear.  As  the  shielding  "charge"  is  built  up  around 
the  tip  with  the  loss  of  the  antishielding  dislocations, 
the  crack  will  require  larger  and  larger  external  A  to 
keep  A,  above  its  critical  cleavage  value.  Ultimately, 
a  steady  state  is  conceivable  in  which  the  shielding 
charge  is  constant,  and  in  this  case,  the  relation 
interconnecting  the  external  A’,  A,*,  w,  m,  etc.,  will  be 
given  by  equation  (27).  Equation  (27)  then  becomes 
a  model  fracture  criterion  for  the  material. 

Equation  (27)  may  be  compared  with  the  elastic- 
enclave  models  of  references  [I]  and  [2],  (We  have 
shown  separately  [13]  that  the  versions  of  the  pre¬ 
vious  elastic  enclave  models  are  equivalent.)  Here, 
however,  we  have  a  more  satisfying  model  of  cracking 
in  terms  of  the  dislocation  theorist.  Our  model,  we 
believe,  should  have  the  same  kind  of  appeal  which 
the  original  BCS  model  did,  and  through  it  we  are 
able  to  see  intuitively  into  the  role  dislocations  and 
their  discrete  sources  play  in  fracture.  We  can  sum¬ 
marize  this  section  by  bringing  together  the  various 
relations  derived  which  interconnect  A,  A,,  it,,  m 

K  =  -2v/2Tt(r,v  R 


o,=  -o„(D*IDar 


D* 


A  I2R 


Our  principal  result,  equation  (27),  is  a  particular 
combination  of  these  equations. 

In  our  model,  the  sources  for  the  dislocations  are 
assumed  to  be  external  to  the  crack  tip,  and  the 
parameter,  tv,  is  a  measure  of  the  inhomogeneity  of 
slip  in  all  real  materials.  If  the  crack  can  itself 
generate  dislocations  out  of  the  tip,  this  process  will 
be  governed  by  an  essentially  different  type  of  physics 
which  we  do  not  address  here. 


CYCLIC  STRESS 

Consider  next  in  a  qualitative  way  "fatigue" 
growth  under  an  applied  cyclic  stress.  Let  the  cyclic 
stress  vary  from  a  =  0  to  a  =  <7ma,  and  the  con¬ 
ventional  stress  intensity  factor  K  vary  from  K  =  0  to 
A  =  Ama,  =  crma, (rt«)' ?.  The  fatigue  crack  will  ad¬ 
vance  an  increment  each  cycle  if  "work  hardening" 
occurs  on  the  two  slip  planes,  both  during  forward 
slip  and  during  reverse  slip.  To  see  why,  consider 
Fig.  7.  Figure  7a  shows  schematically  a  dislocation 
distribution  fl,(.v)  =  flma,(.v)  that  exists  when 
<r  =  (Tmax  and  the  crack  has  stopped  its  advance.  This 
distribution  produces  a  shielding  stress  intensity  fac¬ 
tor  /,ma,  given  by  equation  (9)  to  be 

+  -U';  (29) 
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Fig.  7.  Cyclic  crack  advance.  Dislocation  distribution  on 
one  slip  plane,  (a)  At  maximum  stress  (solid  line)  and  when 
stress  is  reduced  to  zero  (dashed  line),  (b)  Dislocation 
distribution  of  (a)  considered  as  the  original  distribution 
and  a  reversed  slip  distribution  B,{  v).  (c)  After  the  stress  is 
again  increased  to  maximum  value  if  the  crack  tip  were  not 
to  advance.  Forward  slip  can  be  considered  to  occur  by 
creation  of  additional  distribution  ht.  (d)  Crack  advance. 
Dislocation  in  inactive  slip  zone  are  now  "stuck"  in  slip 
plane  and  do  not  advance,  (e)  Activation  of  a  slip  zone  and 
distribution  function  B  containing  both  negative  and  posi¬ 
tive  dislocations.  The  negative  dislocation  eventually  also 
become  stuck,  (f)  After  many  stress  cycles  and  many 
incremental  crack  advances.  A  wake  of  stuck  dislocations  of 
alternating  sign  are  left  behind. 


where  /ma<  is  the  value  of  /  from  (7)  when 
fl(.v)  =  flm,„(.\ ).  Because  the  crack  has  stopped  but 
was  advancing  immediately  before  the  stress  reached 
its  maximum  value 

=  ^nr.i,  +  A,.,;,,  =  A,/,  (30) 

where  K,h  is  given  by  equation  (2b)  and  is  the  critical 
stress  intensity  factor  of  a  Griffith  crack.  Further 
crack  advance  would  reduce  A',  to  a  value  smaller 
than  A', h,  the  critical  value  for  brittle  propagation  of 
the  crack  tip.  Now  let  the  external  stress,  o.  and  stress 
intensity  factor.  A,  be  reduced  to  zero.  Then,  in  the 
absence  of  the  repulsive  forces  on  the  dislocations  due 
to  the  crack  A' -field,  the  large  dislocation  density  at 
the  maximum  near  the  crack  tip  (see  Fig.  6a)  will  be 
forced  in  the  negative  direction  on  the  slip  plane 
(reversed  slip)  under  the  mutual  repulsions  of  the 
dislocations  until  a  new  equilibrium  distribution  is 
achieved,  as  shown  schematically  by  the  dashed  line 
of  Fig.  7a.  If  the  crack  tip  is  locked  in  position,  we 
expect  many  of  these  dislocations  to  be  pushed  well 
behind  the  crack  tip.  When  they  get  into  this  region, 
however,  they  are  in  the  presence  of  the  open  cleavage 
plane  of  the  crack.  Under  these  conditions,  an  inverse 
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image  dislocation  distribution  is  generated  across  the 
open  crack  plane  which  cancels  the  long  range  force 
fields  of  these  relocations.  In  a  sense,  the  crack  faces, 
by  converting  the  whole  dislocations  into  dipoles 
absorbs  these  dislocations,  so  far  as  their  long  range 
stress  fields  are  concerned,  almost  as  effectively  as  if 
they  were  physically  absorbed  into  the  open  cleavage 
faces.  Only  those  dislocations  which  can  be  sustained 
by  the  modified  a,  (see  below)  under  the  mutual 
repulsive  interactions  from  the  dislocations  will  re¬ 
main  in  front  of  the  crack. 

Secondly,  on  a  microscopic  level,  we  can  expect  the 
effect  on  the  crack  tip  of  decreasing  K  to  zero  to  be 
important.  In  the  absence  of  the  external  stress,  <r, 
because  A.',  =  L  from  equation  (30).  K,  is  a  large 
negative  number,  because  L  is  negative.  If  the  crack 
tip  is  reversible,  it  will  thus  be  driven  backwards  on 
its  cleavage  plane.  Of  course,  cracks  in  materials  of 
even  modest  ductility  should  be  reversible  only  in  a 
limited  degree,  and  over  very  short  distances,  but  if 
the  crack  tip  can  recede,  this  provides  additional 
room  for  dislocations  to  stay  ahead  of  the  tip  and  not 
get  lost  behind  it.  In  any  case,  the  negative  K,  caused 
by  the  dislocation  shielding  will  provide  a  large 
closure  stress  on  the  crack  faces,  and  if  the  tip  itself 
remains  open,  the  situation  will  be  analytically 
complex. 

With  this  general  physical  picture  in  mind,  we  turn 
to  an  analysis  of  our  simple  model.  For  definiteness, 
we  assume  the  crack  tip  to  be  locked  in  position,  so 
far  as  retrograde  motion  is  concerned.  After  the 
external  stress  is  turned  off,  but  before  the  dis¬ 
locations  are  allowed  to  redistribute,  K,  =  Z.ma<  <  0,  as 
noted  above.  Because  of  the  reversed  slip,  however, 
the  magnitude  of  L  will  be  reduced.  The  reversed  slip 
can  be  considered  to  occur  by  the  creation  of  a  new 
dislocation  "reversed”  distribution  function  Br(x) 
shown  in  Fig.  7b,  which  when  added  to  the  distribu¬ 
tion  flnm(.v)  gives  the  dashed  distribution  of  Fig.  7a. 
Of  course.  fl,(.v)  must  satisfy  Jflr(.v)d.v  =  0.  The 
value  of  L  is  given  after  reversed  slip  by 


L  —  —  A'na,  +  —  /mjv  —  !,)  ‘ 

(31) 

where 

I  =  1»  MTf  i  vi  it  > 

lU) 

and  Tr(.v)  is  the  total  change  of  stress  that  must  be 
produced  in  the  reversed  slip  zone  at  point  ,v  in  order 
for  slip  to  occur.  For  the  case  in  which  there  is  no 
work  hardening  rr(.v  >  would  be  equal  to  2  <7,„  because 
the  stress  started  out  as  equal  to  a  constant  friction 
stress,  cr0,  before  reversed  slip  started  and  the  stress 
has  to  equal  —  <j„  during  the  reversed  slip.  This  is  the 
value  of  t, (.v ),  however,  only  in  those  regions  of 
reversed  slip  where  0.  Where  flma<(.v)  =  0 

and  reversed  slip  occurs,  t,(.y)  is  smaller  than  2ff„ 
because  the  stress  on  the  slip  plane  was  smaller  than 
<r„  where  flmax(.v)  =  0.  Thus  for  the  case  where  no 
work  hardening  occurs,  r,  is  not  a  constant  over  the 


integral  of  equation  (32).  and  /,  is  not  equal  to  zero. 
Because  B,(x)  is  positive  in  the  regions  (see  Fig  7b) 
where  r,(.v)  is  smaller  than  2  n„  the  term  /,  is  a 
negative  quantity.  The  shielding  intensity  factor  /.  of 
equation  (31)  is  a  less  negative  quantity  (less 
shielding)  than  This  result  is  also  true  when 
work  hardening  occurs. 

Let  the  stress,  which  has  been  reduced  to  zero, 
increase  again  to  the  level  n„u> .  Forward  slip  will 
occur,  as  shown  in  Fig.  7c,  in  the  reversed  slip  zone. 
If  the  crack  tip  were  not  to  advance  the  shielding 
intensity  factor  L*  is  now  given  by 

L  =  -  Amjv  +  (K  ,,uv  —  —  /,  —  //)'■  (33) 

where 

//  =  ^  j\(v  )£,(*)  d.v  (34) 

where  B,(x)  is  the  new  additional  dislocation  distri¬ 
bution  produced  by  the  forward  slip  (JZJ,(.v)djr  =  0) 
and  tffv)  is  the  total  change  in  stress  on  the  slip  plane 
in  the  region  where  B,(x)  ¥=  0.  When  no  work  hard¬ 
ening  occurs  over  the  region  where  flma,(jc)  #  0,  then 
Zf(.x)  =  2  (T0-  But  where  Bmd,(x)  =  0,  the  stress  t,(.y) 
will  be  smaller  than  2  an  by  the  identical  amount  that 
rr(.v)  in  the  same  region  was  smaller  than  2  aa. 
Because  B,{x)=  -  Br(x )  when  no  work  hardening 
occurs,  lf  =  -  /, .  Consequently,  L  =  Lma,  and  no 
crack  advance  will  occur. 

Suppose  work  hardening  occurs  during  slip.  Let 
equation  (25)  be  generalized  for  reversed  slip  (see 
Fig.  5)  to  give  for  the  friction  stress 

ff,  =  (<r„/DJ,)(X|DT,)  +  ff'  (35) 

where  the  sum  on  D*  is  the  sum  of  the  absolute 
values  of  all  displacements  that  have  occurred  at  any 
particular  point  .v  on  a  slip  plane.  (In  reversed  slip  the 
friction  stress  is  given  by  equation  (35),  but  with  a 
negative  sign  on  the  righthand  side  of  the  equation.) 

Because  of  work  hardening  B,(x )  ^  —  B,(x)  and 

^  —Ir  The  increase  in  the  value  of  r,(.v)  and  r,(.v) 
produced  by  work  hardening  can  be  expected  to 
increase  the  magnitude  of  l,  and  t,  given  by  equations 
(32)  and  (34).  However,  the  same  work  hardening  will 
make  the  magnitude  of  Br(x)  and  fl,(.v)  smaller,  an 
.(Iki1  which  wiLUVc'cav  <tu  vzlgevol/  a u«l  /  (Mtv 
power  exponent  m  of  equation  (35)  is  smaller  than  I 
the  effect  of  work  hardening  on  rr  and  r,  will  be 
smaller  than  on  B,  and  B,  and  consequently  \l,  \  <  /,. 
Thus,  if  the  crack  tip  were  not  to  advance  upon 
increasing  the  stress  again  to  a  =  <rmin  the  shielding 
intensity  factor  L  is  smaller  in  magnitude  than  Lma,. 
The  crack  tip  stress  intensity  factor  K,  will  be  larger 
than  Kih.  Thus,  the  crack  tip  must  advance  before  a 
attains  its  maximum  value.  As  the  crack  advances  it 
will  leave  behind  it,  “stuck"  on  the  slip  plane,  part  of 
its  dislocation  distribution  as  indicated  in  Fig.  7d. 
Part  of  the  distribution  is  left  behind  because  on  the 
trailing  edge  of  the  distribution  \B,  \  <  B, . 

As  the  crack  advances  and  spreads  out  the  distribu- 
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tion  J9|(.v)  the  crack  tip  stress  intensity  factor  K,  will 
actually  increase  further  without  increasing  the  stress. 
Note  that  /ma<  is  equal  to 

L.  =~r-  |r(.v)fl,(.v)d.v  (36) 

where  r(.v)  is  equal  to  the  yield  stress  on  the  slip  plane 
wherever  fl,(.v)^0.  Here  fl,(.v)  is  the  distribution 
shown  by  the  solid  line  of  Fig.  7a.  Now  if  the  crack 
advances  and  spreads  out  the  distribution  fl,(.v),  the 
trailing  "stuck"  edge  of  this  distribution  will  be  in 
regions  where  the  stress  on  the  slip  plane  is  smaller 
than  the  yield  stress.  In  the  regions  where  slip  is 
occurring  r(.v)  is  equal  to  the  yield  stress  and  has  the 
same  value  as  the  yield  stress  in  equation  (36).  Thus 
the  "stretched”  dislocation  distribution  BSIr  of  Fig.  7d 
will  lead  to  a  value  of  /slr,  with  sufficient  crack 
advance,  that  is  smaller  than  /m„.  Therefore,  K, 
increases  in  value  with  crack  advance.  But  this  in¬ 
crease  cannot  go  on  indefinitely.  A  large  value  of  K, 
means  that  the  stress  field  of  the  crack  tip  is  increased. 
Eventually  this  increased  stress  field  will  produce  slip 
near,  but  ahead,  of  the  crack  tip  as  shown  in  Fig.  7e. 
This  slip  can  be  described  with  still  another  dis¬ 
location  distribution  function  B'(x)  where  again 
]B'(x)Ax  =0.  The  positive  dislocations  of  this  func¬ 
tion  are  further  from  the  crack  tip  as  shown  in  Fig. 
7e,  and  move  in  the  direction  of  the  advancing  crack 
tip.  The  negative  dislocations  move  in  the  opposite 
direction.  Eventually  the  negative  dislocations  in  turn 
will  also  become  stuck  on  the  slip  plane,  but  positive 
dislocations  will  join  the  moving  dislocations  of  the 
distribution  fl,(.v).  Hence,  the  total  number  of  posi¬ 
tive  dislocations  increases  and  the  value  of  /  also 
increases.  This  increase  in  turn  decreases  the  value  of 
K,  until  once  again  K,  is  equal  to  Kch.  When  the 
applied  stress  a  has  reached  the  value  <Tma> ,  a  limited 
amount  of  crack  growth  thus  will  insure  that  K,  =  Kch 
and  the  fatigue  crack  will  stop.  In  the  wake  of  the 
fatigue  crack  (on  both  of  the  slip  planes)  "stuck" 
dislocations  of  alternating  sign  will  exist  as  shown  in 
Fig.  7f. 

A  rough  estimate  of  the  crack  growth  increment 
per  stress  cycle  can  be  made  as  follows.  To  have  crack 
advance  in  each  cycle,  the  crack  tip  stress  intensity 
factor  K,  given  by  equation  (24)  must  equal  Kth 
during  the  incremental  growth.  The  work  hardening 
under  cyclic  slip  must,  therefore,  increase  the  average 
value  of  the  friction  stress  a,  to  the  point  that 
K,  =  Kch.  Thus  from  equation  (24)  a,  must  attain  the 
value 

ol  =  K,Jf!sfi-  (37) 

where  /(  is  given  by  equation  (27).  Let  fa  be  the 
increment  of  crack  growth  per  cycle.  From  equation 
(20)  when  R  P  w,  the  length  of  a  slip  zone  with  an 
avefage  frieiiofl  stress  o,  is  equal  to  The 

number  of  cycles  n  required  for  the  crack  to  advance 
this  distance  is 


The  average  displacement  D  within  this  slip  zone 
from  equation  (21)  can  be  shown  to  be  (for  one  of 
the  slip  zones  of  the  two  slip  planes)  equal  to 
jf  A'2,a<  6  a,//.  The  sum  1 1 D  |  thus  is  of  the  order  of 
2 nD  or 

I|D|  5  .Jf  JtA'l,  24  a ,’/<fa.  (39) 

Inserting  equation  (39)  into  (35)  gives  for  a, 

(T,  =  (a„jD”)' . . I.jUttA'L,  24 far'* .  (40) 

when  the  term  a'  is  relatively  unimportant. 

Inserting  equation  (40)  into  (37)  gives  the  following 
equation  for  the  incremental  crack  growth  6a  per 
cycle 

fa  =  (Jfjr/fL,  24 //)(<?„/(  v  u  KlhD”fm  * (41) 

This  equation  is  a  fourth  power  Paris  equation. 
Since  the  fatigue  crack  growth  is  essentially  deter¬ 
mined  by  an  accumulated  displacement  criterion,  a 
fourth  power  is  what  is  to  be  expected  [14],  Since  the 
term  Kch  =  (4  /r//Jf)' 2  equation  (41)  predicts  that  any 
environment  that  reduces  the  surface  energy  will 
increase  the  crack  growth  rate. 

It  should  be  emphasized  that  equation  (41)  is  not 
valid  if  either  Kmy  is  so  large  that  fracture  occurs 
under  a  monotonic  increase  of  the  stress  towards  the 
value  Equation  (41)  also  is  not  valid  of  KmdX  is 
so  small  that  the  friction  stress  is  never  reached  on  the 
slip  planes.  If  no  slip  occurs  on  these  planes  the  crack 
will  not  propagate  in  fatigue.  Thus,  if  KmaJ(2nw)' 2  is 
smaller  than  the  friction  stress,  clearly  no  fatigue 
crack  propagation  can  take  place.  In  other  words,  a 
threshold  effect  exists. 


DISCUSSION 

The  original  impetus  for  this  paper  was  to  see  if  a 
difference  between  the  published  results  of  the  first 
and  last  authors  [1,2]  could  be  resolved.  In  one  of 
these  earlier  calculations,  the  fracture  stress  calcu¬ 
lated  for  the  elastic  enclave  model  was  shown  to  be 
proportional  to  y‘ 2  where  y  is  the  surface  energy.  In 
the  other,  the  fracture  stress  contained  in  a  more 
complex  power  relation.  In  further  work,  one  of  us 
[10]  has  also  shown  that  this  discrepancy  is  due  to  the 
way  one  estimates  the  ratio  of  the  plastic  zone  radius 
to  the  elastic  enclave  radius.  However,  in  all  calcu¬ 
lations  of  those  sorts,  in  which  continuum  plasticity 
is  used  in  a  cut-off  procedure  at  the  elastic  zone 
boundary,  very  severe  problems  are  encountered 
regarding  the  kind  of  boundary  condition  to  use  at 
the  elastic  enclave  radius.  In  this  paper,  we  have  thus 
adopted  a  relatively  self-consistent  one  dimensional 
dislocation  model  where  the  boundary  condition  at 
the  "elastic  enclave"  radius  is  not  a  problem. 

Thfi  invKStipltirn  -sport^-H  h  <••-<»  shows  thn»  for  the 
one  dimensional  model,  the  fracture  stress  <rf  is  the 
complex  power  law,  equation  (27) 

K,={K;h!\\)"  1 


n  =nK:mijSa-,Sa. 


(38) 
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Here.  K,  =  o,  no  is  the  critical  stress  intensity  factor 
at  fracture,  and  a  is  the  half  length  of  (he  crack. 
(When  the  Griffith  fracture  condition  determines  the 
value  of  K,h.  this  term  is  equal  to  v4/iy  .*  .)  Un¬ 
fortunately.  however,  this  result  is  valid  only  for  the 
one  dimensional  model,  and  the  full  two  (and  perhaps 
mi vc!  dimensionality  ol  die  silp  uiOu  id  me  ci.uA  up 
might  yield  a  still  different  result. 

One  of  the  features  of  this  work  is  the  modification 
of  the  classical  BC'S  model  of  one  dimensional  slip  in 
front  of  a  crack  to  include  an  elastic  enclave,  or 

disl|>i'.’iti|)n  tfee  /one  Our  calculation  though  inde¬ 
pendently  arrived  at,  is  a  limiting  form  of  other  work 
in  press  by  Majumdar  and  Burns.  This  important 
truth  -tfinWHi  fjtw  ihr  i’JAc'k  lip  tln.'w  imrntity  I'wtii! 
K,  can  he  related  to  the  overall  stress  intensity  factor. 
K.  and  it  allows  us  to  derive  a  fracture  stress  in  terms 
of  materials  parameters  for  the  one  dimensional 
model. 

One  of  the  interesting  side  features  of  this  in¬ 
vestigation  was  the  discovery  that  all  continuum 
theories  such  as  BCS  or  ours  neglects  the  self-image 
term  of  the  dislocations  in  the  crack.  This  self-image 
term,  which  is  present  only  in  bonatide  discrete 
calculations  is  important  in  determining  the  intrinsic 
elastic  enclave  region  around  the  crack  tip,  which 
must  otherwise  be  an  unknown  parameter,  u,  in  the 
theory.  In  our  work,  ic  is  either  a  scale  parameter 
determined  by  the  crystalline  heterogeneity  of  plastic 
flaw  in  the  material,  or  as  a  lower  limit,  it  is  given  by 
the  self-image  term  calculable  in  a  discrete  dislocation 
model. 

In  our  model,  we  have  also  demonstrated  one  of 
the  curiosities  of  the  dislocation  shielding  theory, 
namely  that  under  certain  conditions  it  is  possible  for 
the  shielding  to  be  zero.  In  particular  for  a  group  of 
dislocations  on  certain  slip  planes  in  equilibrium  with 
a  certain  kind  of  friction  stress  of  which  <r,  =  const, 
is  a  special  case,  and  where  the  total  summed  burgers 
vector  is  zero,  the  shielding  is  zero. 

Finally,  we  have  applied  these  ideas  to  cyclic 
fatigue.  In  qualitative  terms  we  showed  that  when  the 
external  stress  cycles  to  “off",  the  large  screening 
dislocation  distribution  near  the  crack  tip  explodes 
under  the  internal  pressure  of  the  mutual  repulsive 
forces  between  the  dislocations.  Many  of  the  dis¬ 
locations  are  pushed  behind  the  crack,  and  "annihil¬ 
ated"  there  by  the  images  in  the  crack  faces,  and  the 
remainder  left  in  front  of  the  crack  is  that  total 
number  which  can  establish  a  new  equilibrium  in  the 


distance  R  under  the  .dip  plane  friction  strew  ir,t  \ ). 
This  residual  shielding  distribution  provides  a  closure 
stress  on  the  lip  which  may  drive  the  crack  tip 
backwards,  and  in  any  case  will  weld  portions  of  the 
cleavage  plane  back  together. 

When  the  external  stress  is  switched  back  "on"  the 
A-iicid  ol  flic  ci etc k  icgc.ie. aics  the  shielding  charge 
partly  by  action  sources  on  the  slip  plane,  and  partly 
by  moving  the  "annihilated"  dislocations  back  into  a 
shielding  configuration  in  front  of  the  crack. 

If  there  is  no  work  hardening,  the  dislocation 
distribution  is  essentially  revervihte  under  (he  stress 
i.e.  the  sources  can  regenerate  exactly  the  same 
distribution  as  in  the  previous  cycle,  with  the  crack  in 
a*  luvf  avtgwi  >.twiii vi  Ihtweirs  if  there  k  -tfrnui 
hardening,  the  sources  are  unable,  at  the  same  stress, 
to  regenerate  the  same  shielding  charge;  K,  rises 
above  Kh.  and  the  crack  moves  forward  to  a  position 
where  fresh  sources  can  be  activated.  Thus  in  our 
view,  fatigue  growth  of  the  crack  is  due  to  strain 
hardening  during  stress  cycling.  Quantitative  analysis 
indicates  that  a  fourth  power  Paris  equation  describes 
the  crack  growth  law. 
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